THIRD ORDER NEWTON'S METHOD FOR ZERNIKE 
POLYNOMIAL ZEROS 



O . RICHARD J. MATHAR 
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' Abstract. The Zernike radial polynomials are a system of orthogonal polyno- 

R ' mials over the unit interval with weight x. They are used as basis functions in 

, ^ ^ optics to expand fields over the cross section of circular pupils. To calculate the 

roots of Zernike polynomials, we optimize the generic iterative numerical New- 
1/^ ' ton's Method that iterates on zeros of functions with third order convergence. 

The technique is based on rewriting the polynomials as Gauss hypergeometric 
functions, reduction of second order derivatives to first order derivatives, and 
. evaluation of some ratios of derivatives by terminating continued fractions. 

' A PARI program and a short table of zeros complete up to polynomials of 

' 20th order are included. 

1. Classical Orthogonal Polynomials: Hofsommer's Newton Method 

^SI ■ The generic third order Newton's Method to compute roots f{x) = numericaUy 

^ \ improves solutions iteratively, starting from initial guesses, 

■ via computation of corrections 

cn 

in 
o 
i> 
o 



- T— I 

>< 



fix)' V 2/'(.t) fix) 

where f{x), f'{x) and f"{x) are the function and its first and second derivatives 
at the current best approximation Xi [H [TTl [M] . For some classes of orthogonal 
polynomials, f" / f can be derived from ///' [l2l I24j. which means the update 
can be done to third order at essentially no additional numerical expense. If we 
divide the differential equation of the classical orthogonal polynomials, for example 



■ as tabulated in [3 22.6] [17], 



(2) h2{x)f" + h^(x)J' + ho{x)f = 

through /', 111) turns into 



(3) A. = -^/ 



1 + + 

2h2(x) f'{x) \ f'{x) J 
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Structure relations 18J relate the ratio ///' to ratios at shifted indices n as tabu- 
lated for example in 22.8], 

(4) g2{x)f'^{x) = gi{x) f n(,x) + gf){x)f n-i{x); 

/gx ^ fnjx) ^ g2{x) 

f'nix) g,{x) + g,{x)^^' 

The benefit is that the three-term recurrence equations, in the notation of ^1, 22.7] 

(6) ai,„_i/„(a;) = (02 ,n — 1 ~r 0-3,71 — ix)fn-l{x) - ai,n-lfn-2{x), 

lead to terminating continued fraction representations for ///' 

fn{x) a2,n-l + a3.n-lX - a4,„ _i ^^'^ 

This is recursively inserted into the denominator of ([5]) to lower the index n until 
/0//1 is reached, which avoids problems with cancellation of digits. 

This work here implements this strategy for the family of Zernikc polynomials, 
/ = i?™, namely (i) fast calculation of /"//' from ///', (ii) calculation of ///' from 
terminating continued fractions, both without evaluation of / or its derivatives via 
direct methods like Horner schemes. 

2. Zernike Polynomials: Derivatives and Roots 

2.1. Definition. We define Zernike radial polynomials in Noll's nomenclature [STJ 
[21[Tll[21[l[2S]forn>0, n-m = mod 2, m < n as 

(ri-m)/2 

(8) Kix) = J2 ~ -~ - 



s=0 



(9) = E (-1) 



s\[{n + m)/2- s\\[{n - to)/2 - s]! 
n — s\ f ri — 2s 



s=0 



s J \{n — m)/2 — s 



Following the original notation, we will not put the upper index m in i?™ — which 
is not a power — into parentheses. The normalization integral is 

(10) f xR^{x)R^,{x)dx= \, J ^,n-. 
Jo 2(n+l) 

The complete set of examples for n < 4 is 

(11) i?°(a;) = 1; RI{x) ^ x; 

(12) i?0(a;) = 2x2-1; RHx) ^ x^; 

(13) Rlix) = 3x^^2x] Rl{x)^x^; 

(14) Rl{x) = 6x^-6x'^ + l; RI{x) ^ Ax^ - 3x^ ; RI{x)^x^. 

The inversion of ^ decomposes powers x^ into sums of i?™(a;), (i > m, i — m even), 

i 

(15) x' = Kn,mRn{x)\ i - m = 0, 2, 4, 6, . . . ; 

n— m mod 2 

(n-m)/2 / -,\s / \ / r, \ 

(—1) I n — s\ I n — 2s \ 



(16) h^^n,m - 2(n+l) J2 



n — 2s + i + 2\ s J\(n — m)/2 — S/ 

S — O ^ / \V // / 
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The basic examples are 










m — 0: 












(17) 


x° = 


Ko[x), 








(18) 


= 


-R'iiix) - 


^ 2 ^\ h 






(19) 


= 


-R'iiix) - 


f -R°^{x) + 


-R°Ax)- 




(20) 


x' = 




^^Rli^) 

20 


+- jRlix) 4 


20 


m — 1: 












(21) 


X — 


ni[x), 








(22) 


x' = 










(23) 


x^ = 


-R](x) ~ 


J_ -Rl(x) + 
t- -"-3 V-^y ^ 






(24) 


x' = 


-R}(x) - 








m = 2: 












(25) 


x" = 


Rl{x)- 








(26) 


x^ = 




f ii?^(a:); 






(27) 


x^ = 


IrH-)- 


f ii?^(.T) + 






(28) 


x' = 


\rI{^)- 


^Y^Rlix) 







Much of this work is based on the representation as a terminating hypergeometric 
function, 

(29) R"^(x) - — F f^H^ _n + rn, J_\ 

with the three negative integer parameters 

(30) a = -{n~m)/2; b = -{n + m)/2; c = -n, 
and the argument 

(31) z^^. 

In these variables, the three-term recurrence reads l9l [23] 

(32) 2n(a-l)(fe-l)i?™+2(x) 

= (n + 1) [2n{n + 2)x^ - rr? - n{n + 2)] i?™(x) - 2ah{n + 2)i?;^'_2(a;). 

i?™(x) is also a product of a;™ times a polynomial of degree n — m, 
(33) 

<'(-) = (-l)^(J'^ + "^J/^)x^ 
\(n — m)i Ij 



a(6- 1) 2 a(a+ l)(b- l)(b- 2) 
m+1 2(m+l)(m + 2) 
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which can be summarized in terms of Jacobi Polynomials 



(34) K{x) 



x"^G-a{m + l,m+\,x'^). 



2.2. Derivatives. Derivatives of (HH) are p!], 3.3.8] 



(35) -K{x) 
ax 

{36)^R^{x) 



(37)^i?-(x) 



n 

{n — m) /2 
n 

{n — m) /2 

dx 
n 

(n - m)/2 



-^x"F{a, b] c; z) + x"-^F(a, 6; c; z) 
dx dx 



-^x^-Fia, b] c; z) + 2-^x"-^F(a, 5; c; z) 
dx- dx dx 



+x"-^F{a,b; c; z) 



d^ d 
— ra;"i^(a, b; c; z) + 3-p^a;"— i^(a, 6; c; z) 
da;'' dx-^ dx 



dx dx'^ 



d^ 



Faa di Bruno's Formula [TUl 0.430.2] [T3] relegates the derivatives w.r.t. x to deriva- 
tives w.r.t. z, 



F{a, b; c; z) 



-j-2F(a,5; c; z) 



i^(a, b; c; z) 



dx 
dx 
dx^ 



-2x'^F'{a,b; c;z); 

6x-'^F'{a, b] c; z) + 4:X'^F"{a, b; c; z); 

-24a;"^F'(a, 6; c; z) - 36a;"^F"(a, b; c; z) - 8x~^F"'(a, 6; c; z). 



After insertion of these three formulas into ([55]) - ([57]) . the derivatives of i?™ = x'^F 
are 

(38) i?"'(2:) = nx''-^F -2x''~^F'; 

(39) i?""(a;) = nin - l)x"-^F ^ 2{2n - 3)x''-*F' + 4x"^^F" ; 

(40) i?""'(a;) = n(n- l)(n-2)a;"-3^-6(n-2)2a;"~^F' 

+12(n - 3)x"-^F" - 8a;"-^i^"', 

where = means the binomial factor and the argument list (a, 6; c; z) of the hypergeo- 
metric function have not been written down explicitly. Since i?™(a;) is a polynomial 
of order n, the (n + l)st derivatives equal zero. Backward elimination of F and its 
derivatives with the aid of [1] 15.5.1] 

(41) z(l - z)F"{a, b;c;z) + [c- {a + b+ l)z]F'(a, 6; c; z) = abF{a, b; c; z) 
leads to the analog of 

(42) .T2(x2-l)^i?™(x)= [nx2(n + 2)-m2]i?™(.T)+x(l-3x2)Ai?™(^), 

which is one special case of differential equations that generate orthogonal functions 
[in], and could also be obtained by applying the derivatives of [2 22.6.1] [SI [3] to 
The derivative of this reaches out to the third derivatives, in which R" is 
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reduced to R and R' with the help of the previous equation, 

(43) x^ix^ - l)2^i?™(x) = [x'^in^ + 2n + 7m^) - 5x''ni2 + n) - 3m^] R'^ix) 

+x [6.t2(2x2 - 1) - m2(a;2 - 1) + 2 + ^^(a;' - l)n{n + 2)] ^Ki^)- 
2.3. Zeros. 

2.3.1. Ratios of Derivatives. Installation of ///' in ^ progresses by dividing _R™ 
x^F through 



(44) 



J?rn'(„\ r, F'(a,b;c;z) 

The analog of ([5]) is implemented by substituting 1, 15.2.1] 

(45) F'(a, b; c; z) = — F(a + 1, 6 + 1; c + 1; z) 

c 

in the denominator. In lieu of ([7|) we find the continued fractions [7] 

(a+l)(c-b)2 (a+2)(c+l-b)z 

^4g^ F(a,b;c]z) _ -bz ^ _^ c(c+i) (c+i)(c+2) 



F(a + 1,5+ l:c+ l:z) c {a+i-b)z -, (a+2-b)z -, 

^ ' ' ' ' c+1 ^ c+2 ^ ^ 

which terminate in our cases since a is a negative integer and c = a + b. This 
already suffices to implement the standard Newton iteration, ie, to approximate 
© by Ax = -fix)/ fix). Division of ^ through i?;"'(a;) yields 



(47) 



Rtix) 



to2\ R^{x) 1 - 3x2 



n(n + 2) 



a;2 J Rl^'ix) 



This is fix)/ fix) of the generic formula, and can be quickly computed from 
i?;j'(a;)/i?™'(x) = fix)/ fix) of the lower order. 

2.3.2. Initial Guesses. For n and m fixed, the strategy adopted here is to compute 
the (n — m)/2 distinct roots in (0, 1) starting with the smallest, then bootstrapping 
the others in naturally increasing order. An approximation to the smallest root is 
found by equating the first three terms in the square bracket of ([33]) with zero — 
hoping that higher powers of x become insignificant for small x — and solving the 
bi-quadratic equation for x. This guess may become unstable for n approximately 
larger than 11 in the sense that the Newton iterations converge to another than 
this smallest root. Instead, the simple, heuristic initial guess 

1.46m + 2.41 

(48) X « 

^ ' n + 0.46m +1.06 

is used in general, but keeping the solution to the bi-quadratic equation when this 
is exact, ie, in the cases n — m = 2 or 4. 

A shooting method is useful to produce an initial estimate for one root supposed 
an adjacent one has already been found. The third order Taylor extrapolation from 
one root x to the next one at x + Ax is 

(49) fix + Ax) « fix) + Ax/'(x) + ^^/" (x) + ^^/'"(x) « 0. 
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Division through Aa;/' and exploiting f{x) ~ yields a quadratic equation for the 
approximate distance Aa; to the next one, 

. , Axf"(x) (Axf f"'(x) 

2 r{x) 6 f'[x) 

from which the branch Ax > is systematically selected to start computation of 
the root adjacent to the previous one. The two ratios of derivatives are obtained 
by setting R^{x) = in (142]) and (|43)). then dividing both equations through 
R^' [x). This aim to locate the next root with sufficient accuracy — and to prevent 
the Newton's Method to be be drawn into the second next root which would call 
for more administrative care [4] — is the rationale to look into third derivatives; it 
might also guide the way to even higher order Newton's methods employing /"'. 

3. Summary 

The Newton's Method of third order convergence is implemented for Zernike 
Polynomials i?™ by computation of the ratios and with relay to 

the generic formulas of associated, terminating hypergeometric functions. Adding 
knowledge on the derivative i?™'", a shooting method is proposed which generates 
an initial guess for the adjacent root from each root found. 

Appendix A. Table of Roots of Low-order Polynomials 

The roots Xi^n,m of are tabulated below for 2 < n < 20 in two major 

columns. Each column contains n, then m, then {n — m)/2 = \a\ values of Xi^n,m- 
Only the roots a; > are included, and only the standard parameter range for even, 
positive values of n — m is considered. 

Because the G'_a(TO+ 1, m+l,y) mentioned in ([M)) build a system of orthogonal 
polynomials with weight over the unit interval Q < y < 1, the squares a;f „ ^ = 
yi,n,m are also the abscissae for Gaussian integration of moment m T, Tab. 25.8] [6]. 



2 





0, 


,7071067811865475727 


14 


4 


0, 


,8246570394661102421 


4 





0, 


, 4597008433809830485 


14 


4 


0, 


, 9278396109654096779 


4 





0, 


,8880738339771152567 


14 


4 


0, 


,9862121569592748882 


6 





0, 


,3357106870197287818 


16 


4 


0, 


, 4352810401596804435 


6 





0, 


,7071067811865475727 


16 


4 


0, 


,6105335011400015999 


6 





0, 


,9419651451198933767 


16 


4 


0, 


, 7533739926853023627 


8 





0, 


, 2634992299855423159 


16 


4 


0, 


, 8653533945767831748 


8 





0, 


, 5744645143153508382 


16 


4 


0, 


, 9446613988337823065 


8 





0, 


, 8185294874300058643 


16 


4 


0, 


, 9894368252983534173 


8 





0, 


, 9646596061808674349 


18 


4 


0, 


,3914606479189805532 


10 





0, 


, 2165873427295972042 


18 


4 


0, 


, 5534677249323020076 


10 





0, 


, 4803804169063914387 


18 


4 


0, 


, 6903188343084598610 


10 





0, 


,7071067811865475727 


18 


4 


0, 


, 8039833969179427386 


10 





0, 


, 8770602345636481223 


18 


4 


0, 


, 8931498697389127495 


10 





0, 


, 9762632447087885579 


18 


4 


0, 


, 9561484720295934103 


12 





0, 


, 1837532119404283737 


18 


4 


0, 


, 9916369448881094950 


12 





0, 


,4115766110542091183 


20 


4 


0, 


, 3554897669109032265 


12 





0, 


,6170011401597257361 


20 


4 


0, 


,5054988151106132310 


12 





0, 


, 7869622564275865484 


20 


4 


0, 


,6353277098425955671 


12 





0, 


,9113751660173390334 


20 


4 


0, 


, 7472059125768437671 


12 





0, 


,9829724091252897145 


20 


4 


0, 


, 8402437401628405356 


14 





0, 


, 1595181614381909196 


20 


4 


0, 


,9130561549207910632 
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14 





0, 


. 3594918736220650279 


20 


4 


0, 


, 9643628549735204780 


14 





0, 


. 5450480935764305812 


20 


4 


0, 


. 9932085369576263423 


14 





0, 


.7071067811865475727 


7 


5 


0, 


.9258200997725514192 


14 





0, 


.8384047803350709316 


9 


5 


0, 


. 7942238940183964369 


14 





0, 


. 9331482158798232174 


9 


5 


0, 


. 9616464847987593600 


14 





0, 


. 9871949939963123866 


11 


5 


0, 


,6840930020506119646 


16 





0, 


.1409080258581174028 


11 


5 


0, 


,8727107021799724862 


16 





0, 


. 3188522562146716699 


11 


5 


0. 


. 9758129824224018867 


16 





0, 


.4870665201405610101 


13 


5 


0, 


.5974058327888663866 


16 





0, 


. 6389700139694938219 


13 


5 


0, 


, 7840085394020771536 


16 





0, 


. 7692316434259740543 


13 


5 


0, 


,9116353458812981314 


16 





0, 


. 8733648750425931917 


13 


5 


0, 


,9831524024611432155 


16 





0, 


. 9478044306220632098 


15 


5 


0. 


. 5288602129223232140 


16 





0, 


.9900226907746954019 


15 


5 


0. 


. 7061939018640647214 


18 





0, 


. 1261740078074202742 


15 


5 


0, 


, 8404366834692439392 


18 





0, 


.2863292621034079777 


15 


5 


0. 


.9345601788393610443 


18 





0, 


. 4396752024502914580 


15 


5 


0. 


.9875190300435001678 


18 





0, 


.5812686885581361818 


17 


5 


0. 


.4737815904024170188 


18 





0, 


.7071067811865475727 


17 


5 


0. 


.6399479279255915198 


18 





0, 


.8137116883158951319 


17 


5 


0, 


,7731681889583568168 


18 





0, 


. 8981568439589463493 


17 


5 


0. 


. 8765920400095927878 


18 





0, 


. 9581312820607194025 


17 


5 


0. 


. 9493949303895553582 


18 





0, 


.9920081248426411147 


17 


5 


0, 


, 9903520628752460198 


20 





0, 


. 1142223084227163565 


19 


5 


0, 


,4287525723551651180 


20 





0, 


. 2597466393536357887 


19 


5 


0, 


,5837872858381362162 


20 





0, 


. 4003688497504367394 


19 


5 


0. 


. 7126994626270685140 


20 





0, 


. 5322614986408245041 


19 


5 


0, 


.8187572351505322255 


20 





0, 


. 6523517689936806363 


19 


5 


0, 


, 9014368896323405878 


20 





0, 


. 7579163340968551044 


19 


5 


0, 


,9596126391703136971 


20 





0, 


. 8465800003925344486 


19 


5 


0, 


,9923041120472186893 


20 





0, 


.9163540713880810040 


8 


6 


0, 


. 9354143466934853324 


20 





0, 


. 9656768006659848247 


10 


6 


0, 


.8164965809277260345 


20 





0, 


. 9934552150241026114 


10 


6 


0, 


, 9660917830792958849 


3 


1 


0, 


. 8164965809277260345 


12 


6 


0, 


. 7132667548676191460 


5 


1 


0, 


. 5958615826865180098 


12 


6 


0, 


. 8854995128634446377 


5 


1 


0, 


. 9192110607898046348 


12 


6 


0, 


.9783359081211411290 


7 


1 


0, 


. 4608042298407784010 


14 


6 


0, 


.6297279581530694781 


7 


1 


0, 


.7684615381131740808 


14 


6 


0, 


.8030111569365681046 


7 


1 


0, 


. 9546790248493448594 


14 


6 


0, 


.9197999176484970008 


9 


1 


0, 


.3738447061866471688 


14 


6 


0. 


. 9847470483288770504 


9 


1 


0, 


.6452980455813291938 


16 


6 


0. 


. 5622633194744699470 


9 


1 


0, 


.8503863747508400017 


16 


6 


0, 


.7291002958196058925 


9 


1 


0, 


.9710282199223060351 


16 


6 


0, 


. 8535798917458220503 


11 


1 


0, 


.3139029878781443572 


16 


6 


0. 


. 9401303969875089983 


11 


1 


0, 


. 5518475574344458012 


16 


6 


0, 


. 9885994756607009437 


11 


1 


0, 


. 7496833930084177977 


18 


6 


0, 


, 5071545248007354179 


11 


1 


0, 


. 8955370355972955831 


18 


6 


0, 


,6650990116151613840 


11 


1 


0, 


. 9798929242261785744 


18 


6 


0, 


,7899637876841181239 


13 


1 


0, 


. 2702856427564344077 


18 


6 


0. 


. 8860816722252334854 


13 


1 


0, 


. 4803812423169180335 


18 


6 


0, 


.9533796195297132847 


13 


1 


0, 


. 6643255837527634045 


18 


6 


0, 


,9911212731459003722 


13 


1 


0, 


. 8142575205172167818 


20 


6 


0, 


,4615059331053839586 


13 


1 


0, 


.9229958831606540626 


20 


6 


0. 


.6100885865453433698 
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13 


1 


0, 


. 9852327505925770890 


20 


6 


0, 


, 7319787736696559133 


15 


1 


0, 


. 2371973029714337655 


20 


6 


0, 


, 8314254884667372503 


15 


1 


0, 


. 4245476318823276363 


20 


6 


0. 


.9085207085711664021 


15 


1 


0, 


.5938221258198196351 


20 


6 


0, 


.9625678871955013483 


15 


1 


0, 


. 7396983346814803850 


20 


6 


0, 


, 9928726539379972849 


15 


1 


0, 


.8568606521572865731 


9 


7 


0, 


. 9428090415820633563 


15 


1 


0, 


.9409149519691435426 


11 


7 


0, 


, 8343946751715023424 


15 


1 


0, 


. 9886964213353295339 


11 


7 


0, 


. 9696128930707756943 


17 


1 


0, 


. 2112674970031354627 


13 


7 


0, 


.7374505105030352281 


17 


1 


0, 


. 3799555987772065824 


13 


7 


0, 


, 8959410734924830866 


17 


1 


0, 


. 5355602273735704522 


13 


7 


0, 


, 9803818873685179458 


17 


1 


0, 


. 674398483981520691 1 


15 


7 


0, 


,6571568797605316092 


17 


1 


0, 


. 7925073093774768207 


15 


7 


0. 


. 8189060737060670503 


17 


1 


0, 


.8863924190810575920 


15 


7 


0. 


.9265772077286310582 


17 


1 


0, 


.9532451899171618948 


15 


7 


0, 


,9860657631628138020 


17 


1 


0, 


.9910701715078688023 


17 


7 


0. 


.5911284719636833440 


19 


1 


0, 


. 1904148441776781775 


17 


7 


0. 


. 7486326613042562084 


19 


1 


0, 


. 3436262195904129513 


17 


7 


0. 


.8647030262108379439 


19 


1 


0, 


. 4870081978929722277 


17 


7 


0, 


. 9448234063345077871 


19 


1 


0, 


. 6179666376955199603 


17 


7 


0, 


,9895076756006223606 


19 


1 


0, 


. 7335050094709572033 


19 


7 


0. 


. 5364117760083304542 


19 


1 


0, 


.8308624555679250401 


19 


7 


0. 


. 6868800285074768697 


19 


1 


0, 


. 9076801990993329516 


19 


7 


0, 


, 8044060078556206639 


19 


1 


0, 


. 9620876312949896425 


19 


7 


0, 


, 8942045500841436789 


19 


1 


0, 


. 9927676416720699892 


19 


7 


0, 


, 9567806932331796022 


4 


2 


0, 


. 8660254037844385966 


19 


7 


0, 


.9917768363310000002 


6 


2 


0, 


.6751652804971347566 


10 


8 


0, 


.9486832980505137680 


6 


2 


0, 


. 9367417879781805290 


12 


8 


0, 


. 8490975736565613552 


8 


2 


0, 


. 5431369901889407936 


12 


8 


0, 


. 9724710674756380513 


8 


2 


0, 


. 8080818238035354373 


14 


8 


0, 


. 7578370705278929531 


8 


2 


0, 


.9628114955311087853 


14 


8 


0, 


. 9046299747290316162 


10 


2 


0, 


. 4518280448392144044 


14 


8 


0, 


.9820745595413366003 


10 


2 


0, 


. 6949479871872660253 


16 


8 


0, 


, 6807458887789931135 


10 


2 


0, 


. 8725819517089140609 


16 


8 


0, 


, 8324053937687126981 


10 


2 


0, 


. 9754483331027854476 


16 


8 


0, 


, 9322943792644380334 


12 


2 


0, 


.3859349518416070879 


16 


8 


0. 


.9871744960883001019 


12 


2 


0, 


. 6047036690554417060 


18 


8 


0. 


.6163473579853351314 


12 


2 


0, 


. 7810976974325815059 


18 


8 


0, 


,7654979183523115127 


12 


2 


0, 


. 9091312772247495122 


18 


8 


0. 


8742424706423237435 


12 


2 


0, 


. 9825584257853499093 


18 


8 


0. 


. 9488319407257733706 


14 


2 


0, 


. 3364437305441 106418 


18 


8 


0. 


. 9902817948430912010 


14 


2 


0, 


.5332155968164153936 


20 


8 


0, 


.5623002252838306125 


14 


2 


0, 


. 7006879382687045688 


20 


8 


0, 


,7059431502741818631 


14 


2 


0, 


.8352024197624841051 


20 


8 


0. 


. 8169643175952393532 


14 


2 


0, 


.9318991674757601817 


20 


8 


0. 


.9012382495289457118 


14 


2 


0, 


. 9869627373443907725 


20 


8 


0, 


, 9597179665033356288 


16 


2 


0, 


. 2980215318345257325 


20 


8 


0, 


, 9923422148430279810 


16 


2 


0, 


.4759020418201234115 


11 


9 


0, 


, 9534625892455923513 


16 


2 


0, 


. 6324385240645979955 


13 


9 


0. 


. 8613939328976448762 


16 


2 


0, 


.7655049675886762550 


13 


9 


0, 


.9748375170636489617 


16 


2 


0, 


. 8714619176716458249 


15 


9 


0, 


,7752632451107936973 


16 


2 


0, 


. 9470538373566472767 


15 


9 


0, 


,9119748762120825081 


16 


2 


0, 


. 9898824305070776930 


15 


9 


0. 


. 9834982084062168228 
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18 


2 


0, 


. 2673780750923551164 


17 


9 


0, 


, 7012604730894272942 


18 


2 


0, 


. 4292182386813452322 


17 


9 


0, 


,8440171326522077910 


18 


2 


0, 


. 5748458299193258680 


17 


9 


0. 


. 9371828916889872740 


18 


2 


0, 


.7031379109502706726 


17 


9 


0, 


.9881197223055366852 


18 


2 


0, 


. 81 13864729725273062 


19 


9 


0, 


,6385852750050469151 


18 


2 


0, 


.8969530821125818454 


19 


9 


0, 


, 7802151464898949840 


18 


2 


0, 


.9576530595755170516 


19 


9 


0, 


,8825163921882515083 


18 


2 


0, 


.9919184646169658670 


19 


9 


0, 


.9522959699168392911 


20 


2 


0, 


. 2423925241972734457 


19 


9 


0, 


.9909494948025596717 


20 


2 


0, 


. 3905933409401172729 


12 


10 


0, 


. 9574271077563381027 


20 


2 


0, 


. 5260526858211680423 


14 


10 


0, 


.8718317153731350855 


20 


2 


0, 


. 6483100675208097741 


14 


10 


0, 


.9768291428674981125 


20 


2 


0, 


.7553451890441794658 


16 


10 


0, 


. 79033474498341 10527 


20 


2 


0, 


.8450496504116072893 


16 


10 


0, 


.9182661541320223941 


20 


2 


0, 


. 9155547407096559231 


16 


10 


0, 


. 9847122732288815516 


20 


2 


0, 


.9653576903957581390 


18 


10 


0, 


. 7192722131418530784 


20 


2 


0, 


.9933952458125041574 


18 


10 


0. 


.8541140545067785750 


5 


3 


0, 


.8944271909999158554 


18 


10 


0. 


.9414112100240797920 


7 


3 


0, 


.7279134123608967943 


18 


10 


0, 


.9889351394546163077 


7 


3 


0, 


. 9480050066727199187 


20 


10 


0, 


. 6583509748272358131 


9 


3 


0, 


. 6027143852742457009 


20 


10 


0, 


.7931750687612736384 


9 


3 


0, 


. 8359493221264154839 


20 


10 


0, 


.8897625152667285597 


9 


3 


0, 


. 9684648164078416555 


20 


10 


0, 


. 9553196436506564693 


11 


3 


0, 


. 5113489892733628084 


20 


10 


0, 


.9915313152497930993 


11 


3 


0, 


.7320153318669290199 


13 


11 


0, 


. 9607689228305228424 


11 


3 


0, 


. 8889787567592866147 


15 


11 


0, 


. 8808037886787085657 


11 


3 


0, 


.9786966233548161087 


15 


11 


0, 


. 9785284968679074380 


13 


3 


0, 


. 4429582456583350258 


17 


11 


0, 


. 8035014494300847243 


13 


3 


0, 


.6458329596901977165 


17 


11 


0, 


.9237159298241287564 


13 


3 


0, 


. 8053384408042754128 


17 


11 


0, 


,9857598779860213822 


13 


3 


0, 


. 9195679148240427647 


19 


11 


0, 


.7352174806998263978 


13 


3 


0, 


.9845992603558400003 


19 


11 


0, 


.8629761939750234534 


15 


3 


0, 


. 3902219391376385849 


19 


11 


0, 


. 9451049233245847336 


15 


3 


0, 


. 5755870542650849409 


19 


11 


0, 


. 9896457777687525104 


15 


3 


0, 


.7296653799715601130 


14 


12 


0, 


.9636241116594315148 


15 


3 


0, 


. 8518232083769392560 


16 


12 


0, 


. 8885993155618584494 


15 


3 


0, 


. 9389444522639546209 


16 


12 


0, 


.9799955389607631906 


15 


3 


0, 


. 9883295764481713208 


18 


12 


0, 


.8151048281585070443 


17 


3 


0, 


. 3484639562887706932 


18 


12 


0, 


. 9284828318578448592 


17 


3 


0, 


.5180380683998889735 


18 


12 


0. 


.9866730715175761057 


17 


3 


0, 


.6639777561205169043 


20 


12 


0. 


. 7494357022184537920 


17 


3 


0, 


.7865138450349050681 


20 


12 


0, 


.8708181789128337513 


17 


3 


0, 


. 8833196205032936010 


20 


12 


0, 


. 9483595986475444883 


17 


3 


0, 


. 9520307422269438380 


20 


12 


0, 


.9902706191583880990 


17 


3 


0, 


.9908430322275690871 


15 


13 


0, 


.9660917830792958849 


19 


3 


0, 


. 3146454476233096487 


17 


13 


0, 


. 8954359308455520639 


19 


3 


0, 


. 4703652913957364068 


17 


13 


0, 


. 9812748728586200286 


19 


3 


0, 


. 6075556764710252633 


19 


13 


0, 


. 8254089266585322715 


19 


3 


0, 


.7270113006012387524 


19 


13 


0, 


. 9326878942873719769 


19 


3 


0, 


. 8270388433632478442 


19 


13 


0, 


.9874761693913842731 


19 


3 


0, 


. 9056952852248679742 


16 


14 


0, 


. 9682458365518542553 


19 


3 


0, 


.9612979313921893310 


18 


14 


0, 


.9014806129319800077 


19 


3 


0, 


. 9926194840347464243 


18 


14 


0. 


.9824003676126956686 
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6 


4 


0, 


.9128709291752769017 


20 


14 


0, 


, 8346212400166618250 


8 


4 


0, 


. 7657261797294159233 


20 


14 


0, 


, 9364250918976404492 


8 


4 


0, 


. 9558574253919850383 


20 


14 


0, 


, 9881879542766879299 


10 


4 


0, 


,6481612911435379321 


17 


15 


0, 


, 9701425001453318764 


10 


4 


0, 


, 8566734238949372804 


19 


15 


0, 


, 9068636367673578169 


10 


4 


0, 


,9726240720110040927 


19 


15 


0, 


, 9833982044520968024 


12 


4 


0, 


,5586908741397313971 


18 


16 


0, 


, 9718253158075500497 


12 


4 


0, 


, 7608650731649617693 


20 


16 


0, 


,9116881483483436632 


12 


4 


0, 


, 9016049670538779370 


20 


16 


0, 


,9842889413994087011 


12 


4 


0, 


, 9811849966449255334 


19 


17 


0, 


, 9733285267845752653 


14 


4 


0, 


,4896856591758441124 


20 


18 


0, 


, 9746794344808964450 


14 


4 


0, 


,6789209633908173114 











Appendix B. PARI implementation 

The full source code of the PARI interpreter program which computed the values 
shown in Appendix \X\ is listed below. The language is similar to C/C++ and has 
inherent support for arbitrary precision computation [35] . 

HypergAugmRatio implements (j46p . HypergRatio implements (j45p . ZernikePrratio 
implements (j44p . Zernike2Prratio implements (I47p . ZernikeSPrratio imple- 
ments (1431) . ZernikeRoot implements ([ij. ZernikeRootEst implements (j48l) . but 
([33|) if n = m + 2. ZernikeAllRoot implements a loop with guesses as in ([50| . 
main loops over n and m to tabulate the zeros up to a maximum n. 

/** Compute the quotient of a Gauss Hypergeometric Function over 

* the Function with the same argument but all the three parameters 

* increased by 1. 

* (§param[in] a first parameter of the Gauss Hypergeometric Function F(a,b;c;z). 

* This must be a negative integer. 

* OparamCin] b second parameter of F 

* (§param[in] c third parameter of F 

* (§param[in] z argument of F 

* Oreturn the ratio F(a,b; c ;z)/F(a+l ,b+l ; c+1 ; z) 

* @warning the function assumes that the parameter a is a negative integer 
*/ 

HypergAugmRatio (a, b, c ,z)={ 
local(ff=0.0) ; 
forstep(ap2= -l,a+l,-l, 

ff = z*(ap2-b)+(ap2+c-a)*(1.0-ff) ; 

ff = ap2*(ap2+c-a-b-l)*z/((ap2+c-a-l)*ff ) ; 

) ; 

return(-b*z/c+l-f f ) ; 

} 

/** Compute the quotient of a Gauss Hypergeometric Function over 

* its first derivative. 

* (§param[in] a first parameter of the Gauss Hypergeometric Function F(a,b;c;z). 

* This must be a negative integer. 

* (§param[in] b second parameter of F 

* (§param[in] c third parameter of F 

* (§param[in] z argument of F 

* @return ratio of the value divided by the first derivative, F/F' . 

* ^warning Checking a against being a negative integer is not done. 
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*/ 

HypergRatio(a,b,c,z)={ 

c*HypergAugmRatio(a,b,c,z)/(a*b) ; 

} 

/** Compute the ratio of the Zernike polynomial over its first derivative. 

* (§param[in] n first parameter of the Zernike polynomial, a positive integer 

* OparamCin] m second parameter of the polynomial, a positive integer 

* less than or equal to n, with n-m even. 

* (§param[in] x the argument of the polynomial in the interval [0,1] 

* @param[in] ffprime the ratio F' (a,b; c;z)/F(a,b; c ;z) of the Gauss 

* hypergeometric function, where a=-(n-m)/2, b=-(n+m)/2, c=a+b, 

* and z=l/x~2. 

* Oreturn the ratio R/R' 

* Owarning no check is done that n-m is a positive even integer 

* or that n aind m are individually positive or that x is between 

* and one. 
*/ 

ZernikePrratio(n,m,x, ffprime) = { 
local(z=1.0/x~2) ; 

/** To enhance stability in the case of ffprime close to 
* zero, we do not use the equivalent x/(n-2*z/ff prime) ; 
*/ 

x*ff prime/ (n*f fprime-2*z) ; 

} 

/** Compute the ratio of the second derivative of the Zernike polynomial 

* over the first derivative, with respect to the argument x. 

* <§param[in] n first parameter of the Zernike polynomial, a positive integer 

* @param[in] m second parameter of the polynomial, a positive integer 

* less than or equal to n, with n-m even. 

* @param[in] x the argument of the polynomial in the interval [0,1] 

* (§param[in] rrprime the ratio R(x)/R'(x) 

* Oreturn the ratio of derivatives R''/R' 

* Owarning no checking is done that n-m is a positive even integer 

* or that n and m are individually positive or that x is between 

* and one. 
*/ 

Zernike2Prratio(n,m,x, rrprime) = { 
local (xsq=x~2) ; 

( (n* (n+2) *xsq-m"2) *rrprime+x* (l-3*xsq) ) / (xsq* (xsq-1) ) ; 

> 

/** Compute the ratio of the third derivative of the Zernike polynomial 

* over the first derivative. Derivatives eire with respect to the airgument x. 

* The ratio R'/R of the first derivative is to be provided. 

* (§param[in] n first parameter of the Zernike polynomial, a positive integer 

* @param[in] m second parameter of the polynomial, a positive integer 

* less than or equal to n, with n-m even. 

* <3param[in] x the argument of the poljmomial in the interval [0,1] 

* (§param[in] rrprime the ratio R'/R 

* Oreturn the ratio R'"/R' 

* Owarning no test is performed that n-m is a positive even integer 
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* or that n and m are individually positive or that x is between 

* and one. 
*/ 

Zernike3Prratio(n,m,x,rrprime) = { 

local (x2=x"2,x2m=x"2-l, n2n=n*(n+2), m2=m"2) ; 
( (x2*(n2n+7*m2)-5*x2"2*n2n-3*m2) *rrpriiiie 

+x* (6*x2* (2*x2-l)-m2*x2m+2+x2*x2m*n2n) ) 
/(x*x2*x2in~2) ; 

} 

/** Compute a root of the Zernike polynomial within a specified error bar. 

* OparamCin] n first parameter of the Zernike polynomial, a positive integer 

* ®parain[in] m second parameter of the polynomial, a positive integer 

* less than or equal to n, with n-m even. 

* @param[in] x an initial guess of the root 

* @param[in] eps the absolute accuracy of the result. The Newton iteration 

* will be terminated if two subsequent estimates agree within this limit. 

* Oreturn a root x such that R(x)=0 

* Owarning no test is performed that n-m is a positive even integer 

* or that n and m are individually positive or that x is between 

* and one . 
*/ 

ZernikeRoot(n,m,x,eps) = -[ 

/** The variable root is used to keep a history of the most recent 

* approximation to the root. The variables a, b, and c are the 

* corresponding arguments to the associated Gauss Hypergeometric Function. 
*/ 

local (root=x,z, a= -(n-m)/2, b= -(n+m)/2, c= -n, fprimef ,ffprime,rrp,r2prp) ; 
if (n==m , 
0.0, 

/** For an absolute accuracy of l.e-30, up to 10 iterations 

* will be needed for n<=10. So we update the root with up to a 

* maximum of 20 Newton iterations. 
*/ 

for(i=l,20, 

z=l/root~2 ; 

/* ffprime contains F' (a,b;c;z)/F(a,b;c;z) of the 

* associated Hypergeometric function. 
*/ 

ffprime= HypergRatio( a, b, c, l/root"2) ; 

/* rrp and r2prp are the ratios R/R' aind R''/R' of 

* the Zernike polynomial relative to its 1st and 2nd derivative. 

* This completes the relay F'/F -> R/R' -> R"/R'. 
*/ 

rrp= ZernikePrratio(n,m,x, ffprime) ; 
r2prp= Zernike2Prratio(n,m,x,rrp) ; 

x=root ; 

/** Now perform the third order Newton update. The formula 

* is Delta(x) = -(f/f ')/[l-(f/f ')*(f "/f ')/2] . Reduction to the 
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* 2nd order would be implemented as the simpler 

* root = x-rrp ; 
*/ 

root = x-rrp/ (1 . 0-0 . 5*rrp*r2prp) ; 

/* terminate if the old and the new guess agree within eps. 
*/ 

if( abs(x-root) < eps, 
break ; 

) ; 

) ; 

retiirn(root) ; 

) ; 

} /* ZemikeRoot */ 

/** Provide a guess of the smallest nonzero root of the Zernike Polynomial. 

* @param[in] n first parameter of the Zernike polynomial, a positive integer 

* @param[in] m second parameter of the polynomial, a positive integer 

* less than or equal to n, with n-m even. 

* @return an estimate of a root which is the smallest nontrivial (nonzero) 

* positive between and 1 out of the total of (n-m)/2 . 
*/ 

ZernikeRootEst (n,m) = { 

local (a= -(n-m)/2, b= -(n+m)/2, x2) ; 
if (n==m, 
0, 

if(n == m+2, 

/** if n equals m+2, the resulting equation is 

* a linear equation in x~2. 
*/ 

x2 = (m+l)/(a*(b-l)) ; 
return ( sqrt(x2)) , 

/* if n does not equal m+2, a simple heuristic estimate 

* with a rational function of n and m is used. 
*/ 

retumC (1 . 46*m+2 . 41) / (n+0 . 46*m+l . 06) ) ; 

) ; 

) ; 

} /* ZernikeRootEst */ 
/++ 

* @param[in] n first parameter of the Zernike polynomial, a positive integer 

* @param[in] m second parameter of the polynomial, a positive integer 

* less than or equal to n, with n-m even. 

* (§param[in] eps the desired absolute accuracy of each root. 

* Oreturn the (n-m)/2 positive roots in the open interval from to 1 

* Onote the m-fold degenerate root at zero is not returned, nor the 

* sjmmietric values on the negative real axis. 

* Qwarning no test is performed that n-m is a positive even integer 

* or that n and m are individually positive or that x is between 

* and one . 
*/ 
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ZernikeAllRoot(n,m,eps) = -[ 

local (s= (n-m) /2 , rs , x , r2pr , r3pr .disc) ; 



/* The result contains s=(n-m)/2 individual values, which are 

* collected in the vector rs. 

*/ 

rs=vector(s) ; 
for(i=l,s, 
if (i==l, 

/** If this is the first root for a pair (n,m) , 

* we attach to the leftmost (positive, but smallest in value) 

* of them, and call ZernikeRootEstO for an initial value. 
*/ 

x=ZernikeRootEst (n,m) , 

/** Otherwise, for the 2nd and higher roots, we taike the 

* previous root, and solve the quadratic equation for the step 

* D to the next root which follows from a 3rd order Taylor 

* approximation at the known root. This estimate is then 

* used to call the Newton routine for the next root . 
*/ 

x=rs[i-l] ; 

/** Compute R''/2R' at the old root. The general formula 

* x~2(x-2-l) [R"/R'] = {n(n+2)x"2-m'2>*[R/R']+x(l-3x-2) 

* is simplified because R and therefore R/R' axe known 

* to be zero for the recent root. 
*/ 

r2pr= Zernike2Prratio(n,m,x,0)/2 ; 

/** Compute R'''/6R' at the old root. Application of the formula 

* x-3(x~2-l)-2 [R' ' VR'] 

* = {x-2(n-2+2n+7m-2. . .>*[R/R']+x*[6x-2(2x-2-l)-m-2(x-2-l)+. . .] 

* uses that R and therefore R/R' are known 

* to be zero for the recent root. 
*/ 

r3pr= Zernike3Prratio(n,m,x,0)/6 ; 

/* The quadratic estimate from l+(r2pr/2)*D+(r3pr/6)*D"2=0 

* is in these variables l+r2pr*D+r3pr*D"2=0, after division 

* l/r3pr + r2pr/r3pr*D +D*2 = 0. disc is the discriminant 

* of the quadratic. 

*/ 

disc = l-4*r3pr/r2pr'2 ; 

/* This third order local approximation to the full polynomial 

* usually has a root to the left and another one to the right of 

* the current root . We select the sign of the square root of 

* the discriminant to lock into the larger of these two values. 
*/ 

if (disOO, 

if( r2pr/r3pr>0, 

X += r2pr/(2*r3pr)*(-l+sqrt(disc)) , 
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X += r2pr/(2*r3pr)*(-l-sqrt(disc)) ; 

) , 

/** If the discriminant was negative, just ignoring the 

* square root works also and keeps arithmetics real-valued. 

*/ 

X -= r2pr/(2*r3pr) ; 

) ; 

) ; 

/** Given the starting value x of the next root, call the Newton 

* routine to converge to the root within eps, and store the result 

* as another value in the vector. 
*/ 

rs[i] = ZernikeRoot(n,m,x,eps) ; 

) ; 

return (rs) ; 
> /* ZernikeAllRoot */ 

/** The main routine tabulates the zeros of the Zernike Polynomials 

* of lowest order. 

* (§param[in] nmax the maximum value of the peirameter n to be used 

* Oreturn 
*/ 

main(ninax)={ 

/** we loop with m and n over all cases that have non-zero roots. 
*/ 

for (m=0,nmax, 

/* The parameter n runs from m up in steps of 2. 

*/ 

f or St ep (n=m , nmax , 2 , 

/* Collect the roots of this (n,m) with an accuracy of l.e-30 

* in the vector rts. 

*/ 

rts=ZernikeAllRoot (n , m , 1 . e-30) ; 

/* Print the (n-m)/2 results. */ 
for(i=l, (n-m)/2, 

print (n" "m" "rts[i]) ; 

) ; 

) ; 

) ; 

return(O) ; 
y /* main */ 

{ 

/** Collect all non-trivial roots up to some maximum explicit n. 

*/ 

main (20) ; 

} 
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